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We have studied the optical oscillation and tunneling of light waves in opti- 
cal waveguide ladders formed by two coupled planar optical waveguide arrays. 
For the band structure, a mid- zone gap is formed due to band hybridization 
and its wavenumber position can be tuned throughout the whole Brillouin 
zone, which is different from the Bragg gap. By imposing a gradient in the 
propagation constant in each array, Bloch-Zener oscillation (BZO) is realized 
with Zener tunneling between the bands occurring at mid zone, which is con- 
trary to the common BZO with tunneling at the center or edge of the Brillouin 
zone. The occurrence of BZO is demonstrated by using the field-evolution anal- 
ysis. The tunable hybridization at mid zone enhances the tunability of BZO 
in the optical waveguide ladders. This work is of general and fundamental 
importance in understanding the coherent phenomena in lattice structures. 
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Light propagation in optical waveguide arrays (OWA) is important both in fundamental 
science and practical applications [1] . Optical waves propagating in a discrete photonic struc- 
ture behave analogously to electrons moving in a semiconductor crystal [1]. Bloch oscillation 
(BO) [2-4] and Zener tunneling (ZT) [5,6] are two fundamental transport phenomena in cur- 
rent research [7,8]. BO is the oscillatory motion of particles in a periodic potential under an 
external force [9], whose optical equivalent is a linear gradient in refractive index [4,8,10,11] 
or geometrical variation in waveguides [12]. In multiband models, ZT acts as the regular out- 
bursts of radiation from BO into higher-order bands [5,13]. The superposition of BO and ZT 
causes a double-periodic motion of wavepackets, which is referred as the Bloch-Zener oscilla- 
tion (BZO) [14,15]. BZO has been predicted theoretically for classical [16] and nonclassical 
light [17] in waveguide arrays and demonstrated experimentally in binary superlattices [8]. 
In most of previous works on BZO, ZT usually occurs at the center or edge of the Brillouin 
zone, that is, tunneling only takes place at specific values of wavenumber. This condition 
restricts the tunability of BZO in applications. To improve the tunability, we analyze the un- 
derlying physical mechanism for the light propagation starting from the dispersion relation, 
because the field evolution is completely controlled by the dispersion relation [18]. 

In this Letter, we consider an optical waveguide ladder (OWL) formed by two coupled 
planar optical waveguide arrays as shown schematically in Fig. 1. The coupling constants 
between neighboring waveguides in the upper and lower arrays are K\ and K2, respectively. 
The inter-array coupling constant is k. In this work, both K\ and k are positive. There 
are two cases for n 2 : (i) k 2 < for the positive/negative coupling (PNC) case and (ii) 
k 2 > for the positive/positive coupling (PPC) case. Although the positive and negative 
coupling constants can be realized by using dielectric and metal OWAs [20-22] , respectively, 
metal OWAs are incompatible with dielectric OWAs due to losses, the PNC is much harder 
to realize experimentally than the PPC, which can be obtained by pure dielectric OWAs. 
However, there is no "true band gap" in the PPC case, while a "true band gap" indeed 
exists in the PNC case. Therefore, we will use the PNC case to illustrate the dynamics of 
BZO, and compare its results with those in PPC case. The numbers of waveguides in the 
upper and lower arrays are the same N up = N\ ow = 100, and the separation between two 
nearest-neighbor waveguides is assumed to be 1. To ensure the occurrence of BO and ZT, 
we impose a linear gradient in the propagation constant, which can be realized by either 
a temperature gradient in thermo-optical waveguide arrays [3, 5] or a gradient in effective 
refractive index by electro-optical effects [2,11]. Light propagates along the waveguide (z) 
axis. 

The propagation of light in the OWL is described by a set of coupled differential equations 
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of modal amplitudes a n (b n ) in the upper (lower) array, 



^dz + ^ an<KZ } + Kl ( a?1 + 1 ^) + a n-i( z )) + K b n (z) = , 
(1-7- + Pi)bn(z) + k 2 [b n +i(z) + b n ^i(z)} + Ka Tl (z) = , 
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where (3 U = (3 + a u n and Pi = Po + 5p + a.\n are the onsite propagation constants of the nth 
waveguide with gradient a u and ai in the upper and lower arrays, respectively We apply 
the same gradient in the two arrays a u = on. The propagation constant difference between 
the two arrays is Sp. The two k terms in Eq. (1) represent the inter-array couplings between 
the two waveguide arrays. Substituting the solutions a™(z) = w™e 1/3mZ and b™(z) = v™e ll3mZ 
into Eq. (1), we obtain a matrix equation H|m) = /3\m), where the Hamiltonian matrix is 
defined as H = {{U, Z}, {Z, V}}. The matrix elements are U i:i = (3 U , U^-x = U iji+ i = k±, 
Vi,i — Pi, Vi.i-i = Vi,i+i = K 2, = k>, and other elements are zero. The column vectors \m) 
and P are eigenvectors and eigenvalues of H, respectively. In the absence of gradient in both 
arrays, a u = a\ = 0, we can obtain a dispersion relation 



The dispersion relation Eq. (2) relates the longitudinal wavenumber P to the transverse 
wavenumber k and determines the evolution of the normal modes during propagation [18]. 
The band structures for PNC case obtained from Eq. (2) are shown in Figs. l(b)-l(d) for 
various values of propagation constant differences Sp between two arrays. The analytic results 
(solid lines) match with the numerical results (dots). There are two minibands /3_ and P + 
with a band gap opening at the mid zone (0 < k < 7r). The physical origin of the mid-zone 
gap is due to the hybridization of the upper array band and the lower array band. The total 
propagation constant difference between the two bands is defined as A = \P + — P_\. The 
band hybridization occurs when the propagation constant difference reaches its minimum, 
that is, dA/dk = 0. The wavenumber k depends on k and 6/3. If there is no coupling between 
two arrays (k = 0), no band gap will be opened. The nonzero inter-array couplings (k ^ 0) 
through evanescent fields between the two arrays can cause the two bands to split at the 
degenerate point due to level anticrossing, which leads to two separated minibands with 
a mid-zone gap. The wavenumber position of the band gap can be tuned by varying the 
propagation constant differences 5p between two arrays. As shown in Figs. 1(b)- 1(d), the 
band hybridization occurs at k = 0.5ir, OAn, and 0.67T when 5p = 0, 1, and —1, respectively. 
In Fig. 1(e), the wavenumber position of the band hybridization covers the whole range of 
[— 7r, 7r] when 5p varies from —3 to 3. For a certain 5p, the wavenumber position of the band 
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hybridization is symmetric with respect to k = 0. The gap due to tunable band hybridization 
is different from the usual band gap appearing at the center or edge of the Brillouin zone 
due to Bragg reflection. The latter is tunable vertically in the gap width, as in binary 
superlattices [8]. While the former can be tuned not only in the gap width but also in the 
wavenumber position. 

In OWL with the obtained band structure, the gap due to band hybridization is so small 
that Zener tunneling between the two minibands is possible. If we impose a gradient in the 
propagation constants in each array, BO and breathing-wave-like oscillation (BW) will occur 
under proper conditions [11]. The superposition of BO and ZT leads to the Bloch-Zener 
oscillation (BZO), and the superposition of BW and ZT leads to the breathing-Zener- wave- 
like oscillation (BZW). The properties of BZO and BZW in OWL will be investigated through 
the field-evolution analysis [11,23]. We inject a Gaussian beam [11] 
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with transverse wavenumber /c and intensity profile \ijj(n, 0)| 2 , which has a discrete Gaus- 
sian distribution centered at the noth waveguide with a spatial width a. The exponen- 
tial factor exp [— iko(n — n )] captures the phase differences between the input beams ex- 
cited at the nth and the n th waveguides. In what follows, we excite the array with a 
plane wavefront with ko = 0. We expand the input Gaussian beam in terms of the su- 
permodes \m) as \i/)(n, 0)) = J2 m Am\ m ), where A m = (m\ip(n,0)) is the expansion am- 
plitude of the input beam. The subsequent wavefunction at a propagation distance z is 
\ip(n,z)) = Ylm A m e ll3mZ \m) . The wavefunction in the reciprocal space can be obtained by 
taking the Fourier transform \(f)(k, z)) = jF\\ij){n, z))\. The evolution of the beam intensity 
\ip(n, z)\ 2 in the real space and \<j)(k,z)\ 2 in the reciprocal space along the propagation dis- 
tance z are used to illustrate the occurrence of BZO and BZW. 

To observe BZO, we use a broad Gaussian beam with spatial width a = 3 to excite the 
OWL from the central waveguide of the lower array. The parameters are a u = ati = a = 0.08, 
6/3 = and k = 0. Figure 2(a) shows the overlapped contour plots of \tf)(n, z)\ 2 in the real 
space for the upper and lower arrays of PNC case, both of which are labeled from 1 to 100. 
Since the coupling constants K\ and in the upper and lower waveguide arrays are positive 
and negative, respectively, the propagation of light beams in the two OWAs is opposite 
in the transverse direction [21]. From Fig. 2(a), the light beam is reconstructed after two 
BO periods, thus the period of BZO Zrzo an d that of BO Z^q satisfy the relationship 
Zbzo = 2Zbo- In the reciprocal space (figures not shown here), ZT occurs at k = 0.5tt and 
k = — 0.57T, which correspond to the spatial positions around Zti = ^bo/4, zt2 = 3Zbo/4, 
zt3 = 5Zbo/4, .... When z < Zti, the light beam mainly undergoes BO in the lower array. 
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At z = zti, ZT occurs, whose tunneling rate is dependent on the width of the band gap. The 
light beam splits into two parts: one tunnels to the upper array, and the other one remains 
in the lower array. Due to the opposite transverse propagation directions in two arrays, the 
two BO paths are different and refocus at the second tunneling point z = z?2- Since both 
bands are occupied, light beams tunnel between each other and coherently interfere. The net 
result is that more field tunnels to the upper array. At z — zt3, more field tunnels back to 
the lower array. After the propagation distance Zz&o, the light beam recovers. Due to the 
occurrence of anomalous BZO, the light beam tunnels at the mid zone between the upper 
and lower arrays periodically. 

It has been proved theoretically [19] and experimentally [5] that light evolution in the 
photonic lattice can reflect the band structure directly, whereas the width of each gap corre- 
sponds to the tunneling rate of ZT [19]. To demonstrate this principle clearly, we take half 
period (Zbo/2) of the light evolution in the real space and rescale by Nh/a for the waveguide 
index and Z B o/(27r) for the propagation distance, as shown in Fig. 2(b). The rescaled light 
evolution matches very well with the band structure described by Eq. (2) shown by the solid 
line in Fig. 1(b) and Fig. 2(b). The BO path with larger oscillation amplitude occurs in 
the upper array, since its coupling constant is stronger which leads to a wider band. ZT is 
demonstrated by the exchange of field intensities of different BO paths at band hybridization, 
where the minimum gap opens. 

For the PPC case, the evolution of broad Gaussian beam in OWL and its visual band are 
shown in Figs. 2(c) and 2(d), respectively. From the band structure of PPC case [solid lines in 
Fig. 2(d)], the two minibands hybridize around k = 0.5tt when 5p = 0, but there is no "true 
band gap" throughout the whole Brillouin zone. Light beam can tunnel between the upper 
and lower arrays around k = ±0.57r, as shown in Fig. 2(c). It indicates that BZO occurs 
in such kind of OWL with pure positive couplings. The different features between the light 
evolution of PPC case [Fig. 2(c)] and PNC case [Fig. 2(a)] are as follows, (i) The transverse 
propagation directions in the two layers of PPC case are the same since their couplings are 
both positive, (ii) The relation between periods of BO and BZO of PPC case is Zbzo — 3Zbo- 
From the visual band for PPC case as shown in Fig. 2(d), we can see that although there 
is no "true band gap", BZO still takes place with tunneling at the wavenumber position 
where the difference between the two minbands at the same k reaches its minimum. This 
phenomena is attributed to the coherent properties of light waves, which is an advantage of 
optical system over other kinds of systems (electronic, plasmonic, elastic, etc). 

For a single waveguide excitation (a = 0.02), as shown in Fig. 3, we observe BZW. 
Figs. 3(a) and 3(b) demonstrate the contour plots of \ip(n, z)\ 2 in the real space for the 
upper and lower arrays, respectively. The light beam enters from the central waveguide of 
the lower array, refocuses at Zbo and recovers at Zbzo — 2Zbo i n the lower array. ZT occurs 
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at Zti = Z B o/4, zt2 = 3Z B o/4, z T3 = 5Z B o/4, ... In the upper array, the light beam cannot 
focus to a single waveguide and can reach a wider range of width due to its stronger coupling 
constant. In comparison, the refocused spatial range of the breathing mode of PPC case is 
much wider in the upper array [Fig. 3(c)] and narrower in the lower array [Fig. 3(d)]. 

The tunabilites of BZO (BZW) are enhanced by the band gap due to tunable hybridiza- 
tion in OWL. Comparing with the band gap in optical binary superlattices [8], which can be 
tuned in width vertically, the band gap in OWL can also be tuned in wavenumber position 
horizontally. Theoretically, we can tune the tunneling position from one array to the other 
one by varying the propagation constant difference 8f3 between two arrays and the tunneling 
rate by the strength of inter-array coupling constant k. The spatial range and intensity of 
field evolution in OWL can also be tuned as well. Moreover, different from the fixed evolution 
pattern in a single OWA, the light beams have different evolution patterns in two arrays in 
OWL and can vary according to different tunneling positions. For the experimental verifica- 
tion, we can resort to the positive/positive coupling case, where the tunable hybridization 
still works and which is much easier to realize through the dielectric OWL. 

In conclusion, we proposed a certain coupled optical waveguide ladder and studied its 
band structure and Bloch-Zener oscillation. Due to the band hybridization through inter- 
array couplings, a band gap opens at the mid zone. This is completely different from the 
conventional Bragg gap. Thus, when a graded potential is applied, BZO (BZW) occurs with 
Zener tunneling at the mid zone. It was found that the tunneling position and rate, spatial 
range and intensity of BZO (BZW) are tunable. These research results are fundamental and 
important in understanding the coherent phenomena in lattice structures. This work will 
arouse interests from two highly active research communities in the fields of lattice systems 
and coherent effects. 
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List of Figure Captions 

Fig. 1. (Color online) (a) Schematic diagram of a cross section of an optical waveguide ladder. 
The coupling constants for the positive/negative coupling case are K\ — 1, K2 — —0.5, and 
k = 0.2. Band structures for various propagation constant differences (b) 8p = 0, (c) 5/3 = 1, 
and (d) 8p = —1. Band hybridization occurs at k = 0.5ir, OAtt, and 0.67T, respectively. 
The analytic results (solid lines) match well with the numerical results (discrete dots), (e) 
Wavenumber position of band hybridization as a function of propagation constant differences 
5p between two arrays. Figl.eps 

Fig. 2. (Color online) Evolution of a broad Gaussian beam centered at ri\ ow = 50 with spatial 
width cr = 3. Positive/negative coupling case (k% = 1.0, k 2 = —0.5, and k = 0.2): (a) Overlap 
of contour plots of \ip(n, z)\ 2 in the n-z domain in the upper and the lower graded waveguide 
array, (b) Visual demonstration of band structure through spatial evolution of light beam. 
Positive/positive coupling case (ki = 1.0, k 2 = 0.2, and k = 0.1): (c) is similar to Fig. 2(a). 
(d) is similar to Fig. 2(b). The common parameters are N up = iVi ow = 100, a = 3, (3 = 2, 
5 p = 0, a = 0.08. Fig2.eps 

Fig. 3. (Color online) Evolution of single waveguide excitation starting from the lower array 
center ni ow = 50. Positive/negative coupling case (k\ = 1.0, k, 2 = —0.5, and k = 0.2): 

(a) Contour plots of \ip(n,z)\ 2 in the n-z domain in the upper graded waveguide array and 

(b) that in the lower graded waveguide array, respectively. Positive/positive coupling case 
(ki = 1.0, K 2 = 0.2, and k = 0.1): (c) is similar to Fig. 3(a), (d) is similar to Fig. 3(b). 
Parameters are the same as those in Fig. 2 except for a = 0.02. Fig3.eps 
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Fig. 1. (Color online) (a) Schematic diagram of a cross section of an optical 
waveguide ladder. The coupling constants for the positive /negative coupling 
case are K\ = 1, «2 = —0.5, and k = 0.2. Band structures for various prop- 
agation constant differences (b) 5p = 0, (c) 5p = 1, and (d) 5p = —1. Band 
hybridization occurs at k = 0.5tt, OAtt, and 0.67T, respectively. The analytic 
results (solid lines) match well with the numerical results (discrete dots), (e) 
Wavenumber position of band hybridization as a function of propagation con- 
stant differences 5p between two arrays. 
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Fig. 2. (Color online) Evolution of a broad Gaussian beam centered at n\ ow = 
50 with spatial width a = 3. Positive/negative coupling case (ki = 1.0, k 2 = 
—0.5, and k = 0.2): (a) Overlap of contour plots of \i/)(n, z)\ 2 in the n-z domain 
in the upper and the lower graded waveguide array, (b) Visual demonstration 
of band structure through spatial evolution of light beam. Positive/positive 
coupling case (ki = 1.0, k 2 = 0.2, and k = 0.1): (c) is similar to Fig. 2(a). (d) 
is similar to Fig. 2(b). The common parameters are N np = Ni ow = 100, a — 3, 
p Q = 2, 5 P = 0, a = 0.08. 
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Fig. 3. (Color online) Evolution of single waveguide excitation starting from 
the lower array center rii ow = 50. Positive/negative coupling case (ki = 1.0, 
K2 = —0.5, and K = 0.2): (a) Contour plots of \ip(n, z)\ 2 in the n-z domain in 
the upper graded waveguide array and (b) that in the lower graded waveguide 
array, respectively. Positive/positive coupling case (ki = 1.0, k 2 = 0.2, and 
k — 0.1): (c) is similar to Fig. 3(a), (d) is similar to Fig. 3(b). Parameters are 
the same as those in Fig. 2 except for a = 0.02. 
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